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ON THE STABILITY OF THREE-DIMENSIONAL BOUNDARY LAYERS 
PART 2: SECONDARY INSTABILITY * 

ERIK JANKE AND PONNAMPALAM BALAKUMAR 

Abstract. The secondary instability of three-dimensional incompressible boundary layers is studied 
using Floquet theory. Starting from the equilibrium solutions that we obtained from the PSE computations 
documented in Part 1, we investigate the region where a purely stationary crossflow disturbance saturates for 
its secondary instability characteristics utilizing developed global and local eigenvalue solvers that are based 
on the Implicitly Restarted Arnoldi Method, and a Newton-Raphson technique, respectively. The main 
focuses of this study are on the existence of multiple roots in the eigenvalue spectrum that could explain 
experimental observations of time- dependent occurrences of an explosive growth of traveling disturbances, 
on the routes by which high-frequency disturbances enter the boundary layer, as well as on gaining more 
information about threshold amplitudes for the growth of secondary disturbances. 

Key words, crossflow, Floquet theory, secondary instability, implicitely restarted Arnoldi method 

Subject classification. Fluid Dynamics 

1. Introduction. In the present work, secondary instability is defined as the linear instability of a 
periodic secondary flow. Due to the nonlinear interaction and saturation of primary stationary disturbances, 
highly inflectional velocity profiles are formed in all coordinate directions as documented in the first part of 
the present paper. According to Orszag and Patera [16], these inflectional profiles are the origin of this in- 
viscid instability to high-frequency disturbances that is inherently three-dimensional. Once the disturbances 
reached a certain threshold amplitude, their explosive growth rapidly leads to transition. 

Herbert [7] presented results from an investigation of the secondary instability in plane channel flow. 
Incorporating the shape- assumption and starting from equilibrium solutions, he solved the linearized dis- 
turbance equations as an eigenvalue problem using Floquet theory. He found that the three-dimensional 
secondary disturbances travel at slightly different phase speeds than the primary two-dimensional distur- 
bances. He further determined threshold amplitudes for the onset of different secondary instability modes 
for primary disturbance amplitudes of less than one percent of the freestream velocity and argued that the 
threshold amplitude physically represents a minimal vorticity concentration that must overcome the viscous 
damping of the secondary instabilities. Finally, he pointed out that there is an intricate connection between 
linear and secondary instability modes that is revealed by considering the limit of a zero amplitude of the 
primary disturbance. 

Reed [14] investigated the existence of a parametric resonance in the region of a finite- amplitude crossflow 
vortex for flow conditions adapted to the swept flat plate experiment by Saric and Yeates [15]. Using the 
parallel flow assumption for the mean flow, neglecting the mean flow distortion as well as the amplitude 
variation in the streamwise direction due to the primary disturbances, and finally superimposing two three- 
dimensional secondary waves in a moving frame of reference on this flow, she established a Floquet system 
that was solved using a spectral collocation method. An unstable structure corresponding to the second 
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harmonic of the primary disturbance was found away from the wall, and the results closely reproduced the 
experimentally observed modification of the dominant wavelength by Saric and Yeates. 

Prom his classical experiment on a swept cylinder, Poll [17] reported occasional instances where high- 
frequency disturbances were measured on top of the primary disturbances. The frequency of these distur- 
bances was an order-of- magnitude higher than the frequency of the primary disturbances. In the experimental 
work by Kohama et al. [9] and by Takagi and Itoh [18], it was pointed out that the high-frequency secondary 
instability originates away from the wall in regions of high vorticity and shear at normal locations of about 
half the boundary layer thickness. Moreover, Deyhle and Bippes [4] mentioned observations of an explosive 
growth of traveling modes at the very end of the transition process that looked like a bursting of individual 
stationary vortices. They also measured a high-frequency disturbance of /*= 2.1 kHz riding on top of the 
primary waves prior to the local breakdown to turbulence. 

Numerical results by Balachandar et al. [1] for rotating-disc flow and by Malik et al. [13] for Swept 
Hiemenz Flow showed that there are several unstable secondary eigenvalues present in the region of nonlinear 
saturation, and that the developing instability structures travel on the back of the primary stationary 
crossflow vortices. Their growth was attributed to combined vortex stretching and tilting in the presence 
of strong flow gradients. Both for rotating-disc flow as for Swept Hiemenz Flow, the inclination of the 
secondary instability structures with respect to the constant phase lines was determined between 40°-50°. 
In contrast to the small threshold amplitudes found by Herbert [7] for two-dimensional boundary layers, 
Balachandar et al. [1] determined threshold amplitudes of about 10% for the onset of secondary instabilities 
in three-dimensional boundary layers. Also in contrast to the transition in two-dimensional boundary layers, 
Balachandar et al. reported similar characteristics for the fundamental and the subharmonic secondary 
disturbance type. 

Wintergerste and Kleiser [19] presented results from a temporal DNS computation of the transition 
process in the DLR Experiment. In the region of nonlinear saturation, they reported the development of 
a secondary crossflow vortex close to the wall, which is in qualitative agreement with the PSE results by 
Malik et al. [13]. Continuing their simulation further into the transition process, Wintergerste and Kleiser 
found that the secondary vortex initiated a vortex-splitting into smaller wavelength structures that caused 
the turbulent breakdown shortly thereafter. 

A qualitatively new result indicating the existence of algebraically amplified modes in the otherwise 
convectively unstable Blasius boundary layer was presented by Koch [8]. Computing nonlinear equilibrium 
solutions at finite Reynolds numbers using numerical bifurcation theory, he examined the secondary insta- 
bility of the such obtained modified mean flow using Floquet theory. In contrast to previous work, the 
nonlinear computation of the modified mean flow allowed him to omit the shape assumption. Investigating 
secondary disturbances that travel with the same phase speed as the primary disturbances (phase- locked), 
Koch found the existence of links between several unstable modes (coalescence) and of a modal degeneracy. 
The latter is defined as the coincidence of two or more eigenvalues and their eigenvectors that leads to a 
locally algebraic growth of these disturbances. It is this mechanism that might explain the experimentally 
observed explosive growth of high-frequency disturbances in the late stages of the transition process. 

In summary of the presented review, the secondary instabilities are seen as high-frequency disturbances 
that originate from regions of high shear and vorticity away from the wall and that grow explosively within 
a short streamwise distance. Their frequency is usually in the kHz-range, which is an order-of- magnitude 
higher than the frequency of the most amplified primary traveling disturbances and in accordance with the 
typical local time scale t* = S*/U* in the region of the nonlinear saturation. The explosive growth of 
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the high-frequency disturbances then initiates a break-up of the large-scale crossflow vortex structures into 
smaller wavelength vortices that is immediately followed by the breakdown to turbulence. From the above, 
it is clear how the secondary instability to high-frequency disturbances develops. Important open questions 
regard threshold amplitudes necessary for a self-sustained growth of the high-frequency disturbances, and the 
route by which the high-frequency disturbances enter the boundary layer before their amplitude is explosively 
amplified by the strong growth rates typical for the secondary instability. 

In the present second part of this paper, we investigate the secondary instability of a mean flow that is 
modified by the presence of a stationary crossflow vortex using Floquet theory. The equilibrium solutions 
for the considered problems of Swept Hiemenz Flow and the DLR Transition Experiment are computed 
using the Parabolized Stability Equations (PSE), as described in Part 1. The formulation of the Floquet 
theory is introduced in Section two. Starting with a brief summary of the observations in Part 1, Section 
three presents results from a temporal analysis of the secondary instability of Swept Hiemenz Flow. Here, 
the emphasis is on the routes by which the secondary disturbances enter the boundary layer, and on their 
threshold amplitudes. Section four contains similar results for the DLR Experiment. There, however, the 
focus is on the existence of multiple eigenvalues in the region of nonlinear amplitude saturation. 

2. Floquet Theory. A Floquet system is defined to be a system of linear ODE’s whose coefficients are 
periodic in the independent variables. Floquet Theory was, for example, applied by Herbert [7] to explain 
the nonlinear transition process in two-dimensional boundary layers, by Reed [14] for a swept flat plate 
boundary layer, as well as by Fischer and Dallmann [5] for the DLR Transition Experiment. 



Fig. 2.1. Coordinate systems and wave angle definitions 


Using Floquet theory, we investigate the stability of a new mean flow consisting of the original basic flow 
and the saturated primary disturbance wave. We define a new Galilean coordinate system (xi,X 2 ,xs) that 
travels with the phase velocity in the xj -direction. The coordinate Xg is aligned with the constant phase lines 
as shown in Figure 2.1. In this coordinate system, the new mean flow Q 2 = {Us, lUg, Us, P&} T is formed 
by a superposition of the parallel base flow Qo(%s) and the solutions of the primary disturbance quantities 
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obtained from a PSE computation. In equation 2.1, a 3 is the real wave number in the x } -direction. Its defi- 

OO 

( 21 ) Qs(xi,x s )=Qo(x s )+ £ g/ m (*s)e im “**' 

m= — oo 


( 2 - 2 ) a s = sgn(<P)Jaf^Tp? ; >P j = tan J ( — ) 

,rea/ 

nition is given in equation 2.2, where is the wave angle of the primary disturbances. By forming the 
mean flow Q 2 as described above, the neglect of the mean flow distortion (shape assumption) and the 
assumption of an arbitrary amplitude of the primary disturbance utilized in previous work ((7), [5]) are 
avoided. The temporal stability of this mean flow is investigated by seeking a disturbance of the form shown 
in equation 2.3. There, the new vector of unknowns q 3 consists of the four primitive variables u 3 , w 3 , v 3 
and p 3 . Further, the real part of cr 3 measures the growth of the secondary instability, its imaginary part 
determines the frequency, and fi 3 is the real wave number tangential to the constant phase lines. 


(2.3) 


OO 

qs(x l ,x s ,x s ,t) = e a ' t+i ' ) ' £ ‘ ^ qs m (x 3 )e ,ma ' £ ‘ 

m= — oo 


Substituting equations 2.1 and 2.3 for V# and vx , respectively, into the Navier-Stokes equations in 
disturbance form (see Part 1), dropping the nonlinear term ( V] * V) Vi and introducing a new vector of 
unknowns A m ={u Srn , dus m /dx 3 , w 3m , dw 3m /dx 3 , v 3m1 Ps m } T , we obtain a linear system of ODE’s for the 
vector of the secondary disturbance eigenfunctions A. The resulting equations for all the Fourier modes are 
given in equations 2.4 - 2.17, where A is a 6 (2*M-f 1) vector at the normal location k . 


(2.4) 
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(2.15) 


mj nj mi rit ^ 

(2.16) A= -^(m 2 a 2 Sm + 0$) + a 3 ; 

(2.17) TT =AA 

Collecting terms with the same Fourier component, writing equations 2.4 - 2.16 in the entire wall- 
normal domain, and employing a fourth-order- accurate compact scheme formulation, we obtain a generalized 
eigenvalue problem of the form given in equation 2.18. In fact, the two-point, fourth-order accurate compact 
scheme formulation applied here ([12]) yields the non-trivial case of a singular, non- symmetric, and semi- 
positive-definite generalized eigenvalue problem. The leading dimensions of the coefficient matrices A and B 
in equation 2.18 can easily reach significant orders. For example, considering a minimal problem size for a 
satisfactory resolution in both the Fourier space and the wall-normal direction, a leading matrix dimension 
of 9-6-71=3834 is obtained by distributing 71 points in the normal direction and truncating the Fourier series 
at M= 4. Hence, the need for efficient eigenvalue solvers is obvious. 

(2.18) AA = o s BA 

Utilizing the QZ-algorithm that is implemented in the ZGEGV-routine available in the public-domain 
software library LA PACK proved to be extremely CPU-time intensive for these matrix dimensions. For 
example, it takes about five CPU-hours on a Sun-Ultra-2 workstation (333 MHz) to compute the entire 
eigenvalue spectrum of a problem with a leading matrix dimension of 2700. Therefore, we adapted the 
recently developed Implicitly Restarted Arnoldi Method ([11]) to the present problem. The method is 
available as a part of the public- domain software library A R PACK. It approximates the eigenvalues and 
eigenvectors in specified regions of the eigenvalue spectrum. Applying this method to a problem size of 5000, 
the required CPU- time to compute 10 eigenvalues that are located in a selected region of interest amounts 
to only one minute on the workstation quoted above. 

For a confirmation of the eigenvalues obtained from the global solver, we developed a local eigenvalue 
solver based on a Newton-Raphson technique. The solution method is based on a two-point, fourth-order- 
accurate compact scheme formulation ([12]), and the solution algorithm allows for several iteration options 
in order to account for the presence of (2-Af+l) eigenmodes in the coefficient matrices. The iteration is 
performed at the wall-normal location where the most dominant secondary eigenmode has its maximum. 
There, we drop either one of the momentum equations, or the continuity equation, solve the bio ck-tri diagonal 
system, and check the convergence by evaluating the dropped momentum, or continuity equation. Hence, 
the developed solution method is capable of iterating on different eigenmodes of the unknowns U 3 , W 3 and 
V 3 at a specified wall-normal location. 

3. Swept Hiemenz Flow. Provided that a mean flow exists that is modified by a saturated crossflow 
vortex, the secondary instability originates from the locations of high shear and vorticity away from the 
wall. This is seen in Figure 3.1 where we plot the vorticity of the total flow quantity Q 2 that is computed 
according to equation 3.1. 

(3.1) \ n \ = 

At the streamwise location of ite=643.5, a high vorticity concentration is seen at normal locations 
between X 3 ~ 1 and Xs ~ 2, which corresponds to about half the boundary layer thickness. Also seen are 
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the footprint of the stationary primary crossflow vortex at the wall that corresponds to the streaks of oil or 
china-clay observed in experiments, as well as a developed secondary crossflow vortex close to the wall that 
shows as A-shape structure at Re= 6 43.5. 



Fig. 3.1. Total vorticity in the presence of a stationary crossflow vortex 

The temporal investigation of the secondary instability is started at the streamwise location where 
#e=546. At this location, the stationary disturbance of an initial amplitude A s = 0.1% has fully saturated and 
reached an amplitude level of A s = 24% for the uj -component (see Figure 3.2 in Part 1). The computations 
presented here are limited to a maximal number of eight modes ( m = -8,* ■ -,8) in the Fourier approximation. 
However, Malik et al. [13] showed that including 16 modes does not change the results qualitatively. Further, 
the results presented consider the fundamental type of secondary disturbances only. This restriction of the 
work is supported by the findings of Balachandar et al. [1] for the three-dimensional boundary layer of a 
rotating-disc flow where they did not find a qualitative difference between the behavior of fundamental and 
subharmonic secondary disturbances. 

The computational grid applied in the present computations is described as follows. For the same 
computational domain as in Part 1, where the freestream boundary conditions were enforced at a normal 
extension of about 10 boundary layer thicknesses, the total number of points is decreased from 141 to 81. 
Also, the grid is now stretched in the entire domain according to equation 3.2 such that the number of points 
within the boundary layer is 50. For the local solution method, the iteration can be performed on either the 
continuity equation, or one of the momentum equations. Additionally, the disturbance mode to be iterated 


(3.2) 


x s = 


A • 77 


1 + 


Xs, rr 


- 7 ? 


on can be specified. Table 1 presents the findings from a grid refinement study. Shown are results from 
the local method for different grids and iteration options. The guess values from the second-order-accurate 
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global method for the two different numbers of points are given by 0 -$== (0.0267 ,-1.2043) for 71 points, and 
<75=(0. 0262, -1.2073) for 81 points. The results of the fourth-order-accurate local method vary only in the 
fourth decimal place, and thus, the grid resolution is considered satisfactory. 


Table 3.1 

Grid study at Re=546, as =-0.508, fl$=0.9, Mode=-l 


A 

Points 

Equation 

3*3, max 

°3 

3 

71 

continuity 

44 

(0.021652,-1.212607) 

3 

71 

x^-mom. 

44 

(0.021652,-1.212607) 

5 

71 

continuity 

44 

(0.021779,-1.212359) 

5 

71 

x^-mom. 

44 

(0.021779,-1.212359) 

3 

81 

continuity 

55 

(0.021715,-1.212264) 

5 

81 

continuity 

55 

(0.021729,-1.212671) 


In general, there will be several unstable eigenmodes for a given wave number at a streamwise location. 
Thus, a scan of the complex cr^-plane for a wave number of (3 3— 0.8 is performed first. Here, the global method 
provides the guess values that are checked for using the local solver. All eigenvalues that were found within 
the scanned domain are plotted in Figure 3.2. Even though they appear clustered around < 7 s = (0.,-l.), there 



Fig. 3.2. Temporal eigenvalue spectrum at Re=546 

is no obvious connection between the positive (unstable) real parts. For the chosen parameters, this isolated 
appearance indicates that there are no multiple eigenvalues present in the flow. The five detected unstable 
eigenvalues are given in the order of increasing frequency in Table 3.2. It is seen that the unstable secondary 
frequencies are an order- of- magnitude higher than the frequency considered in Part 1 (F=0.7510 -4 ). 

The shape- functions of the five unstable modes are shown in Figure 3.3. Plotted is the summation of 
eight Fourier components of the wave number a . 3 in the Galilean coordinate system. It is observed that 
the shape-functions of the different modes vary widely in features like the location of the maximum in the 
normal direction, the largest disturbance component, the magnitude, the presence of multiple peaks, and the 
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Table 3.2 

Unstable eigenvalues Re=546 , as =-0.508, fls=0.8 


Mode 

os 

F — O 8 ,imag / Rz 

1 

(0.0113,-0.5160) 

1.03-10 -3 

2 

(0.0030,-0.8566) 

1.71-10 -3 

3 

(0.0467,-0.9595) 

1.9210- 3 

4 

(0.0218,-1.0838) 

2.1610- 3 

5 

(0.0021,-1.1205) 

2.2410 -3 


extension into the outer flow. Of special interest are Modes 4 and 5, as well as the most unstable Mode 3. 
Because of their full shape- profiles, Modes 4 and 5 might be easier detected in experiments than the other 
three modes that show a more narrow and fluctuating shape. For Mode 4, in particular, it is noted that the 
component tangential to the constant phase fines has a shape that is very similar to the experimentally and 
numerically observed profile in the DLR Transition experiment at xj/c=0.80 (Fischer and Dallmann [5]). 
For that reason, most of the following studies consider Mode 4. 

In Figure 3,4, the shape- functions of the individual Fourier components are given for Mode 4. Plotted are 
the normalized quantities for the first seven modes. It is observed that the disturbance component tangential 
to the constant phase fines (0-Mode) shows a double peak structure for both wall-parallel disturbance 
quantities. The maximal disturbance quantity is found for the component at x$ — 3, where its normalized 
Value \^S } o\/\^S ,total\max — 1 • 

While studying the origin of the secondary instability, other important questions address the existence of 
a fink between the eigenmodes of the undisturbed mean flow, investigated by the linear stability theory, and 
the unstable secondary eigenmodes, investigated by Floquet theory, and also the magnitude of the threshold 
amplitudes at which the secondary disturbances become unstable. In order to investigate these questions, 
the amplitude of the stationary disturbances obtained from the nonlinear stability analysis in the previous 
section is set to ,4=1.0 for the given Reynolds number. Then, by gradually decreasing A to zero, the variation 
of the eigenvalues is monitored. During the computations, the evolution of the shape- functions is followed in 
order to assure that no jumps to different eigenvalues occur. A typical step size in the amplitude is ,4A=0.01, 
and the computations are continued until the local method ceases to converge. 

For a Reynolds number of ite=546, a value of A=1 corresponds to amplitudes of A s = 24%, A a = 17% and 
A 9 =l% for the u*-, and vi -disturbance components of the primary stationary vortex, respectively. In 
Figure 3.5, the variation of the eigenvalues crs-ios^ai, &3,imag) versus a decreasing amplitude A is shown 
for the five unstable eigenvalues at A= 1 . Additionally, the eigenvalues of the undisturbed parallel mean 
flow Qo(xs) are shown in the relevant <7 5 -range. It turns out that only one of the five unstable eigenvalues 
at i?e=546 has a fink to the eigenvalue spectrum of the undisturbed mean flow. This connection could be 
established for Mode 1, the mode with the lowest unstable secondary frequency (F=1.0310 -3 ). As seen in 
Figure 3.5, the variation of the eigenvalues with decreasing amplitude A takes very different routes. For 
example, the highest frequency mode (Mode 5) becomes stable at a threshold amplitude of A thrcs =0.98. 
Shortly thereafter, the iteration for an eigenvalue ceases to converge. On the other hand, Modes 2 and 3 
could be followed further until A~0.1. The threshold amplitudes for Modes 1-4, respectively, are A ^^=0.88, 
AtAres=0.95, A ^*=0.31, Athres—' 0-55. Thus, this analysis predicts the onset of a secondary instability due 
to Mode 3 already for stationary disturbance amplitudes of A*=0.31-24%=7.5% and A s — 0.31*17%=5% for 
the Ui- and Wi -components, respectively. The corresponding values for Mode 4 are A a =0.55-24%=13% and 


8 




\/\ws\ max \ws\l\w3\ max M/|u>s| max 


Fig. 3.3. Total shape-functions for unstable modes at Re=546 , as =‘0.508, fis =0.8 


A s = 0.55*17%=9%. 

In order to investigate the eigenmode variation in the amplitude range where the eigenvalue computa- 
tion stops converging, the next study examines the decay of the individual disturbance components at the 
freestream boundary. Thereby, it can be determined whether the convergence problems are due to an insuf- 
ficient grid resolution, or if they are part of the physics of the investigated problem. In the present stability 
formulation, a complex eigenvalue reaches the continuous spectrum when the corresponding eigenfunction 
does not decay in the freestream and is of rather oscillatory nature ([2]). In the following, the method of the 
investigation is described, and results are presented for the Modes 1 and 4. 

In Floquet theory, we seek a solution to the linear eigenvalue problem that is written as in equation 2.17 
(see Section two), where the elements of the coefficient matrix A are determined from a nonlinear PSE 
computation. Since the coefficients of the matrix are not constant in the entire domain, one usually solves 
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Fig. 3.4. Shapes of the secondary instability eigenfunctions for Mode 4 at Re=546 


the eigenvalue problem using a Newton-Raphson technique. Here, however, the focus is on the eigenvalues of 
A at the freestream boundary where the primary disturbance quantities are zero, and thus, the coefficients of 
the matrix A are constant. Therefore, the behavior of the secondary disturbances at the freestream boundary 
can be investigated by finding the eigenvalues of the submatrix A n at the maximal normal extension of 
the computational domain. The dimension of the submatrix A n depends on the number of considered 
Fourier modes in the Floquet analysis and is determined to (2«Af+l)-6 for the six components of the vector 
of unknowns (see Section two). Considering eight Fourier modes in the analysis, the complex and non- 
symmetric coefficient matrix A n consists of 102x102 elements. If the eigenvalue problem to be solved is 
well-posed, one obtains 51 positive and 51 negative eigenvalues which are distributed almost symmetrically 
about the imaginary axis. In Figure 3.6a), the eigenvalue spectrum at the freestream boundary of Mode 4 
is shown for an amplitude of 4=1. In that case, none of the eigenvalues has a zero real part, the secondary 
disturbances decay in the freestream, and the iteration for an eigenvalue of the entire problem converges. 


Table 3.3 

Eigenvalue characteristics at the freestream boundary for Mode 4 


A 

1.0 

0.8 

0.6 

0.4 

0.38 

0.36 

0.34 

#0/Areoi,i > 0 

51 

51 

51 

51 

51 

51 

53 

$ Of ^real,i < 0 

51 

51 

51 

51 

51 

51 

49 

(^reoi,*)mtn 

0.8 

0.8 

0.8 

-0.16 

-0.078 

-0.011 

0.0066 


When approaching the continuous spectrum, however, the real part of the eigenvalue closest to the 
imaginary axis tends towards zero, and the well-posedness of the eigenvalue problem is violated by an 
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Fig. 3.5. Growth rates and frequencies of the secondary disturbances vs. the amplitude A (Re=546, as — -0.508, {3$ =0.8 ) 


a) 


b) 




Fig. 3.6. Eigenvalue spectra at the freestream boundary for Mode 4 ( a) A=l, b) A =0.34) 


imbalance of the eigenvalues with a positive and negative part. This is seen in Figure 3.6b), where the 
eigenvalue spectrum at the freestream boundary of Mode 4 is shown for an amplitude of ;4:^0.34. There 
are two eigenvalues close to the imaginary axis which indicates that the disturbances do not decay in the 
freestream. Further, tables 3 and 4 list the number of eigenvalues with positive and negative real parts in 
the spectrum, as well as the magnitude of the eigenvalues closest to the imaginary axis for the Modes 1 and 
4 in dependence of the amplitude A. It is seen that the eigenvalues of Mode 4 join the continuous spectrum 
at an amplitude of ^4=0.34, whereas the Mode 1 indeed reaches the linear eigenvalue spectrum, as stated 
earlier. 

In Figure 3.7, the normalized total shape- functions of the ^-component of Mode 1 are given for different 
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Table 3.4 

Eigenvalue characteristics at the freestream boundary for Mode 1 


A 

1.0 

0.8 

0.6 

0.4 

0.2 

0.1 

0.01 

#of\ real , x > 0 

51 

51 

51 

51 

51 

51 

51 

ifcof ^rea/,t ^ 0 

51 

51 

51 

51 

51 

51 

51 

(^real,i)min 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 


amplitudes. It is seen that the shape-functions for the different amplitudes indeed belong to the same family 
of eigenmodes. Comparing the Unear shape-function with the nonlinear shape-function at >1=0.01 , it is 
obvious that nonlinear and linear results merge in the limit of A — ► (). 



a) A=1.00 b) >1=0.90 c) >1=0.80 



d) >1=0.70 


e) >1=0.60 


f) >1=0.50 



g) .4—0.25 h) .4=0.01 i) Linear 


Fig. 3.7. Shape-Junctions |ti) 3 |/|tZi 9 )mox | of the linear eigenmode and of Mode 1 for various amplitudes A 


Next, the eigenvalues (r s for the different unstable modes are computed as functions of the wave number 
tangential to the stationary vortex at Re= 546. In particular, three of the five unstable eigenmodes at 0 S = 0.8 
are traced through their unstable wave number range. Figures 3.8 and 3.9 show the growth rate and frequency 
development, respectively, that was obtained by closely watching the evolution of the eigenfunctions with 
changing wave numbers. Thereby, it is assured that the curves for the different modes in Figure 3.8 indeed 
belong to the specified modes. 

From Figures 3.8 and 3.9, the following observations can be made. A wide range of unstable secondary 
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Fig. 3.8. Secondary growth rates at Re =546 



Fig. 3.9. Secondary frequencies at Re=546 

wave numbers is present at the investigated Reynolds number. Starting at /?$=(), which corresponds to an 
alignment with the primary wave vector, the disturbances turn towards the direction of the constant phase 
lines with a growing wave number tangential to the stationary vortex. The dependence of the frequency on 
the wave number /3s is approximately linear, which corresponds to a constant phase speed of the disturbances. 
An intricate structure of the high- frequency /low- growth rate Mode 5 is observed. The growth rate behaves 
in a periodic manner as a fraction of the wave number (3$ . The most unstable disturbances of Modes 4 and 
5 are found in the wave number range of $5 =0.6 - 0.9. This corresponds to an inclination of the secondary 
structures of ^=-30° -40° with respect to the direction of the constant phase lines. The Modes 4 

and 5 have very similar frequencies; however, the frequency curves never cross each other. Even though 
the curves for the real parts (growth rates) of the different modes do cross each other, their eigenvalues 
Gs—iGStveaU &s,irnag) do not coincide because of the missing corresponding crossovers of the imaginary parts 
(frequencies). Thus, multiple eigenvalues were not detected at the investigated location. 
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In order to gain more insight into the spatial development of the secondary instability, a temporal 
analysis at eight different streamwise positions is performed. Choosing Mode 4 at lle=546 and decreasing 
the Reynolds number, this mode is traced back to its onset at Re^ 475. There, the amplitude level of both 
the uj- and the ui ; -components of the stationary vortex is .4, =11%, which closely corresponds to the results 
obtained by decreasing the amplitude A for Mode 4 (see Figure 3.5). Figure 3.10 shows the secondary growth 
rates at the different streamwise positions. It can be seen that the unstable spanwise wave numbers align 
more closely with the direction of the wave vector for a decreasing Reynolds number. 



Fig. 3.10. Secondary growth rates at different streamwise stations (Mode 4) 



Fig. 3.11. Trace-back of the secondary instability (Mode 4) 

In Figure 3.11, the most amplified secondary growth rates are plotted together with the primary growth 
rates. The most unstable secondary frequencies are found in the range of F=1.010~ 3 -2.010 -3 , which is 
an order-of-magnitude higher than the most amplified primary frequency of F=1.0-10~ 4 . Also, a linear 
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development of the temporal growth rate crs, real is observed. Table 5 shows the wave numbers a 3 , the 
secondary growth rates, the wave numbers 0s , the wave angles of the secondary disturbances according to 
the definition in Section two, as well as the dimensional frequencies of the points of maximal amplification 
at the considered Reynolds number locations. 


Table 3.5 

Parameters for the most unstable disturbances of Mode 4 


Re 

as 

® 3 ,real,max 

03 ,as=max 

^ 2 ,cts —max 

F„ 3 = max ■ 10 3 

546.0 

-0.508 

0.0218 

0.80 

-32.4° 

2.167 

538.5 

-0.510 

0.0192 

0.70 

-36.1° 

1.886 

523.5 

-0.513 

0.0142 

0.50 

-45.7° 

1.332 

516.0 

-0.515 

0.0123 

0.50 

-45.8° 

1.305 

508.5 

-0.517 

0.0102 

0.40 

-52.3° 

1.028 

501.0 

-0.519 

0.0081 

0.45 

-49.1° 

1.144 

486.0 

-0.524 

0.0030 

0.40 

-52.6° 

0.987 

478.5 

-0.527 

0.0006 

0.40 

-52.8° 

0.973 


Another interesting feature in Figure 3.11 is the development of the primary and secondary growth 
rates. At ite~510, the primary and secondary growth rates are of the same order. At i?e=546, however, the 
primary growth rates are about to become stable, whereas the temporal theory predicts a continued steep 
growth of the high-frequency disturbances. Looking at Figure 3.1 in Part 1 where the amplitude evolution of 
the stationary disturbances was plotted, we note that the onset of the secondary instability according to the 
temporal theory occurs at a location where the stationary disturbances start to saturate, and both the uj- 
and the wi -components of the stationary primary disturbance have reached an amplitude level of A s = 11%. 

4. DLR Transition Experiment. In this section, results from a temporal analysis using Floquet 
theory performed at the chordwise location of £j/c= 0.6 are presented. At this location, the averaged 
stationary disturbance component tangential to the inviscid streamline Utang.avg with an initial amplitude 
of A s ~0.1% at xijc— 0.1 has reached an amplitude level of 20% and is almost saturated. Special 
attention will be directed towards the existence of multiple eigenvalues in the considered region of nonlinear 
saturation that might explain the experimentally observed time-dependent occurrences of an explosive growth 
of traveling disturbances ([4]). 

Pursuing the sequence of global and local eigenvalue computations as described in Section three, we find 
that there are more unstable eigenvalues present at this location than detected for Swept Hiemenz Flow 
at the investigated locations. The results for the secondary growth rates and frequencies of three selected 
modes are plotted in Figures 4.1 and 4.2. The unstable eigenvalues appear densely clustered in the wave 
number range between /3?=0.1 and /?$=0.4. The disturbances with these wave numbers have frequencies 
of /*=500 Hz up to /*=4000 Hz (see Figure 4.2). The development of the eigenfunctions with a varying 
wave number tangential to the constant phase lines 0$ was followed closely in order to assure that the 
plotted eigenvalues indeed belong to the same family of unstable secondary eigenmodes. In contrast to 
Swept Hiemenz Flow, where neither multiple eigenvalues, nor crossing growth rate curves could be found at 
the investigated locations, it can be seen in Figure 4.1 that the growth rate curves of the different modes 
cross. This indicates that the algebraically unstable modes found by Koch [8] for Blasius boundary layer 
flow might be also present in this three-dimensional boundary layer flow. Table 6 lists the secondary growth 
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rates, the wave numbers /3 S , the wave angles of the secondary disturbances according to the definition in 
Section two, as well as the dimensional frequencies at the points of maximal amplification for the considered 
modes. It is observed that the wave angle range of the most amplified secondary disturbances agrees with 
the results for Swept Hiemenz Flow. 



00 0.2 0.4 0.6 0.8 

P 3 


Fig. 4.1. Secondary growth rates at xi/c=0.6 



Fig. 4.2. Secondary frequencies at xj /c=0.6 

Further, Figure 4.3 shows the total shape- functions of the most unstable secondary modes at xj /c=0.6. 
Plotted are the normalized shape- functions of the disturbance velocity components for the most amplified 
wave numbers f3 $ . Comparing the shape- functions of the uns table modes in the Swept Hiemenz Flow problem 
at 546 and /3$=0.8 (see Figure 3.3) with the shape- functions shown in Figure 4.3, the following is noted. 
First, the shape- functions at Re = 546 are fuller than the profiles obtained here. This is attributed to the 
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Table 4.1 

Parameters for the most unstable modes at xijc=0.6 , (as =-0.201) 



& S , real , max 

$3 s~max 

^2 ,<j $ = max 

r 1 

J( 7 S —max ! 

Mode 1 

0.018001 

0.35 

-29.9° 

3117 

Mode 2 

0.032616 

0.35 

-29.9° 

2933 

Mode 3 

0.001906 

0.19 

-46.6° 

1451 


higher amplitude level of the stationary disturbances in the Swept Hiemenz Flow problem (A 3>u =24% at 
Re= 546; A 3 jU =20% at X\jc~ 0.6). Second, all the shape- functions at Xijc— 0.6 show a strong maximum at 
a height of x$ — 1 mm, which is equivalent to about 30% of the boundary layer thickness at this location 
and can be interpreted as a formation of shear layers near X 3 ~ 1 mm. 



\u 3 \/\w s \ max | | / 1 Wg | rnax \j \^3 \ max 

Fig. 4.3. Total shape-functions of the most unstable modes at xj /c=0.6, (as =-0,201, normal axis in millimeters) 


From Figure 4.2 and as observed for Swept Hiemenz Flow, it appears that a linear frequency dependence 
on the wavenumber can serve as a first guess for a distinction of the different modes. For closely clustered 
eigensolutions, however, only an additional comparison of the eigenfunctions can be a conclusive test. It 
is also noted from Figure 4.2 that the frequency curves of Modes 1 and 3 cross at the same wave number 
#9=0.15 as the corresponding growth rate curves in Figure 4.1 (see Point 1 .). Hence, the existence of 
a multiple eigenvalue at this wave number needs further investigation. The temporal wave number at the 
location where Modes 1 and 3 cross is < 75 =(0.00146,-0.33792). This corresponds to a dimensional frequency of 
/*=1145 Hz, which is about six times the most amplified frequency. In contrast, the dimensional frequencies 
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at the <Ts -locations where Modes 1 and 2 cross at /?$= 0.5 are/*=4561 Hz and /*=4291 Hz for Mode 1 and 2, 
respectively. Even though the growth rates are similar at this wave number, the difference in the frequency 
is considered too large to suspect a multiple root there. 

In Figures 4.4 - 4.6, the eigenfunctions of the two crossing Modes 1 and 3 of Figure 4.1 are plotted 
at three wave numbers at and near the point of their crossover. It is seen that the shape- functions of all 
disturbance quantities are very similar close to the crossover location at /? 5 =0.15. Further, the development 
of the shape- functions of both modes with an increasing wave number tangential to the constant phase lines 
indicates that the two modes indeed coincide. Hence, the existence of a multiple eigenvalue at xi/c=0.6 can 
be concluded. 






Fig. 4.4. Total shape- functions near the multiple eigenvalue at xj /c=0.6 ((3 S =0.125) 






Fig. 4.5. Total shape-functions near the multiple eigenvalue at Xi fc=0.6 (/3s =0.150) 

Concluding the secondary instability analysis of the DLR Transition experiment, the structure of the 
secondary instability is examined. Considering the mode with the highest frequency at xj/c= 0.60 (Mode 1), 
Figure 4.7 shows contours of the amplitude of the secondary eigenfunction tangential to the constant phase 
lines |ws| for the most amplified wave number /3 5 =0.35, and iso-lines of the modified mean flow component 
in that direction in the Galilean coordinate system. 

Superimposing the secondary structure on a modified mean flow component Ws that is assumed to be 
constant in the xj -direction, the total flow component tangential to the constant phase fines Wg — W s + w 5 
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Fig. 4.6. Total shape-functions near the multiple eigenvalue at xj/c=0.6 ((3s —0.190) 



Fig. 4.7. Contours of |u;$| and W 2 at xifc=0.6 (axis in mm) 


is obtained. Shown in Figure 4.9 is a top view of Ws contours at the same wall-normal location of x<?= 2.2 as 
in Figure 4.8, where the secondary structure is assigned an amplitude of 1%. The dark patches correspond 
to lower-speed fluid traveling on the crest of the primary crossflow vortices in the x^-direction. The phase 
speed of these secondary disturbances C p hase = \&s,imag\/ Ps * U^q/U^ is C phase — 1- Depending on their 
assigned amplitude, these secondary structures can severely modify the ordered primary structure of the 
stationary crossflow vortex. 

From the observations in Figures 4.7 - 4.9, it seems that the structure of a single- frequency secondary 
instability is well defined and might be detected in experiments, if care is taken as to where the measurements 
are performed in the boundary layer. 

5. Conclusion. In the second part of the present paper, we studied the secondary instability originating 
from a three-dimensional mean flow that is modified by the presence of purely stationary disturbances. 
Starting from equilibrium solutions that were obtained from nonlinear PSE computations, we used Floquet 


19 








Fig. 4.8. Contours of w s at xj/c=0.6 and x s =2.2mm (as =-0.201, Ps=0.S5) 


theory and developed global and local solvers that are based on the Implicitly Restarted Arnoldi Method 
and a Newton-Raphson technique, respectively. 

For both the Swept Hiemenz Flow problem and the DLR Transition Experiment, several unstable eigen- 
values were detected in the considered region of nonlinear amplitude saturation. The frequency of these 
traveling secondary disturbances is an order-of-magmtude higher than for the most amplified traveling pri- 
mary disturbances. In contrast to the Swept Hiemenz Flow problem, where the few detected unstable and 
isolated eigenvalues indicate the dominance of a convective instability type in the transition process, the 
existence of a multiple eigenvalue was documented for the DLR Experiment at a chordwise location of 
Xx/c= 0 . 60 . The observation of this locally algebraic instability in an otherwise convectively unstable flow 
might explain the experimentally observed phenomena of a time dependent explosive growth of traveling 
disturbances. 

For Swept Hiemenz Flow, the existence of a fink between the unstable secondary eigenvalues and both 
the eigenvalue spectrum of the undisturbed mean flow and the continuous spectrum was established. It is 
concluded that there are at least three different mechanisms in the region of nonlinear saturation that cause 
high-frequency secondary instabilities. First, due to the action of the growing nonlinear disturbances, stable 
linear eigenmodes of the undisturbed mean flow are modified such that they develop into unstable secondary 
eigenmodes. Second, an inviscid high-frequency instability originates from the highly inflectional character of 
the modified mean flow profiles. Third, due to a not yet investigated receptivity mechanism, high-frequency 
disturbances from the continuous spectrum present in the outer flow enter the strongly distorted boundary 
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Fig. 4.9. Contours of Ws at xj fc=0.6 and xs— 2.2mm 


layer in the region of nonlinear amplitude saturation. This corresponds to observations by Choudhari et 
al. [3], where they allocated the receptivity mechanism to regions of a strong variation in the boundary 
layer thickness. However, it is understood that the second mechanism represents the prevailing path to the 
transition in a crossflow instability dominated boundary layer. 

Further considering Swept Hiemenz Flow, an unstable secondary mode at Re = 546 was traced back to its 
critical point by gradually decreasing the Reynolds number. The onset of this secondary instability was found 
at Re ^ 475 in a region where the primary stationary disturbances are about to saturate. The amplitudes 
of both the ui - and Wj -disturbance components of the stationary crossflow vortex at this Reynolds number 
are A s = 11%. This value agrees well with the prediction of the threshold amplitudes necessary for a self- 
sustained growth of the secondary disturbances of Mode 4 as obtained by artificially decreasing the primary 
disturbance amplitudes at Re= 546 ( A s = 13%). Other detected eigenvalues at Re = 546 become unstable at 
threshold amplitudes of the primary disturbances as low as A s =5% (Mode 3). Hence, secondary instabilities 
with small growth rates can be assumed to be present in a crossflow instability dominated flow even before 
the primary disturbances start to saturate. The impact of these weak secondary instabilities on the overall 
transition process, however, needs further investigation. 

For both investigated flows, the obtained wave angles of the most unstable secondary disturbances 
closely correspond to previous work ([1], [13]). The most amplified secondary structures are inclined at 
angles between ^=-30° and #£=-55° with respect to the direction of the constant phase lines. Future 
work related to the research presented in Parts 1 and 2 will particularly focus on modeling the secondary 
instability more straightforwardly using the nonlinear PSE. 
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